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ANALYSIS  AND  COMPARISON  WITH  THEORY  OF  FLOW-FIELD 
MEASUREMENTS  NEAR  A  LUTING  ROTOR  IN 
THE  LANGLEY  FULL-SCAIE  TUNNEL 
By  Harry  H.  Hey son 

SUMMARY 


Results  of  an  investigation  in  the  Langley  full-scale  tunnel  of  the 
induced  flow  near  a  lifting  rotor  are  given.  Measurements  of  stream 
angles  and  velocities  were  made  in  several  transverse  planes  along  and 
behind  the  rotor  in  four  different  conditions  representative  of  the 
cruising  and  high-speed  ranges  of  flight.  These  measurements  indicate 
that  available  theory  may  be  used  to  calculate  with  reasonable  accuracy 
the  induced  flow  over  the  forward  three-quarters  of  the  disk  for  these 
flight  conditions  provided  that  a  realistic  nonuniform  rotor  disk-load 
distribution  is  assumed.  Rearward  of  the  three-quarter-diameter  point, 
calculations  of  the  induced  velocity  are  increasingly  inaccurate  due  to 
the  rolling  up  of  the  trailing-vortex  system.  Farther  rearward,  well 
behind  the  rotor,  the  flow  may  be  represented  more  accurately  by  the 
flow  behind  a  uniformly  loaded  wing. 


INTRODUCTION 


A  detailed  knowledge  of  the  induced  velocity  in  and  near  the  rotor 
is  required  for  a  complete  helicopter  analysis.  Fortunately,  very  super¬ 
ficial  estimates  are  usually  adequate  for  isolated  rotor  performance  or 
stability  calculations.  Recently,  however,  the  emphasis  on  rotor-blade 
vibration  problems  and  the  advent  of  the  compound  helicopter  and  converti- 
plane  have  highlighted  the  need  for  more  complete  information. 

Rotor-blade  vibration  analyses,  in  general,  require  a  knowledge  of 
the  instantaneous  induced  velocity  acting  on  the  rotor  blades  at  any 
particular  moment.  Some  experimental  evidence  (ref.  l)  is  available 
which  indicates  that  the  fluctuating  component  of  the  flow  may  be  so 
large  that  there  is  little  apparent  relation  between  the  instantaneous 
and  the  time-averaged  induced  velocities.  Such  instantaneous  velocities 
have  not  been  calculated  theoretically  and  are  not  discussed  in  this 
paper . 
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Although  performance  calculations  for  a  compound  helicopter  or 
convertiplane  do  not  require  as  detailed  a  knowledge  of  the  flow  as  do 
vibration  analyses,  they  do  require  a  more  elaborate  analysis  than  is 
needed  for  an  isolated  rotor.  A  knowledge  of  the  time-averaged  induced 
velocity  at  locations  in  the  vicinity  of  the  rotor  should  be  adequate  to 
predict  the  effect  of  the  rotor  on  the  wings  and  tail  surfaces  of  these 
machines.  The  time-averaged  flows  are  considerably  simpler  to  deal  with 
theoretically  than  are  the  instantaneous  flows,  and  the  results  of  several 
theoretical  analyse-?*  are  available  (refs.  2  to  6,  for  example).  These 
analyses  are  limited  by  the  restricted  number  of  locations  for  which  the 
calculations  were  made  and  also  by  the  fact  that  the  side -to- side  asym¬ 
metries  caused  by  a  finite  tip-speed  ratio  are  usually  neglected. 

Several  previous  experimental  investigations  (see  ref.  7  for  a 
bibliography)  have  provided  some  qualitative  information  on  the  subject 
by  the  use  of  smoke  in  flight  and  in  wind  tunnels.  However,  quantitative 
measurements  (refs.  8  to  10)  are  extremely  scarce  and,  in  general,  are 
inadequate  either  to  define  the  nature  of  the  flow  or  to  check  the 
accuracy  of  the  theoretical  calculations. 

This  investigation  was  undertaken  in  an  attempt  to  provide  such 
information.  Stream  angle  and  velocity  measurements  were  made  in  several 
transverse  planes  along  and  behind  a  lifting  rotor  in  the  Langley  full- 
scale  tunnel.  The  simulated  flight  conditions  cover  tip-speed  ratios 
between  0.095  and  0.232  and  calculated  wake  skew  angles  between  75-0° 
and  85.8°.  Substantially  all  the  basic  data  are  presented  in  order  to 
provide  the  designer  with  the  maximum  amount  of  information.  At  those 
locations  for  which  theoretical  calculations  are  available,  the  measured 
data  are  compared  with  theory  in  order  to  determine  the  extent  to  which 
the  calculations  are  valid  and  also  to  determine  the  general  nature  of 
the  flow  field. 


SYMBOLS 


rotor  lift  coefficient, 


C>p  rotor  thrust  coefficient, 


p(fiR)2nR2 


equivalent  flat-plate  area  representing  parasite  drag,  based  on 

unit  drag  coefficient.  Helicopter  parasite  tog  n 

PV^ 
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i  complement  of  skew  angle  (used  in  ref.  4  as  angle  of  attack  of 

rotor  tip-path  plane),  deg 

L  lift,  lb 

q  dynamic  pressure,  lb/sq  ft 

( 

free-stream  dynamic  pressure,  lb/sq  ft 

R  blade  radius,  ft 

s  semi span  of  wing,  ft 

T  rotor  thrust,  lb 

V  forward  speed  of  rotor,  ft/sec 

v  normal  component  of  induced  velocity,  positive  downward,  ft/sec 

vq  average  or  momentum  value  of  the  normal  component  of  induced 
velocity,  positive  downward,  ft/sec 

X  distance  parallel  to  longitudinal  rotor  tip -path-plane  axis, 

measured  positive  rearward  from  center  of  rotor,  ft  (fig.  l) 

Y  distance  parallel  to  lateral  rotor  tip-path -plane  axis,  measured 

positive  on  advancing  side  of  disk  from  center  of  rotor,  ft 
(fig.  1) 

Z  distance  parallel  to  axis  perpendicular  to  rotor  tip-path  plane, 

measured  positive  above  rotor  from  center  of  rotor,  ft  (fig.  l) 

a  rotor  control-axis  angle  of  attack,  radians 

T  circulation,  sq  ft/sec 

7  vorticity,  radians/sec 

M-  tip-speed  ratio,  — 

uR 

p  mass  density  of  air,  slugs/cu  ft 

X  rotor  wake  skew  angle,  angle  between  Z-axis  of  tip-path  plane 
and  axis  of  skewed  wake,  positive  rearward  from  Z-axis,  deg 
(fig.  1) 

n  rotor  angular  velocity,  radians/sec 
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Note  that  the  definitions  of  X,  Y,  and  Z,  together  with  the  direc¬ 
tion  of  rotation  of  the  rotor  in  this  investigation,  result  in  a  left- 
handed  axis  system. 


APPARATUS  AND  TESTS 

The  surveys  were  conducted  in  the  Langley  full-scale  tunnel  which 
is  described  in  reference  11.  The  experimental  apparatus  used  in  these 
tests  is  shown  in  figure  2(a). 

The  rotor  was  of  the  teetering  or  seesaw  type  and  had  untapered, 
untwisted  blades  (fig.  2(b))  with  NACA  0012  airfoil  section.  The  rotor 
radius  was  7*5  feet  and  its  solidity  was  0.0543.  The  rotor  tip  speed 
was  500  feet  per  second  for  all  tests  except  that  at  |i  =  0.232.  For 
this  test,  the  rotor  was  operated  at  450  feet  per  second. 

Two  different  survey  rakes  were  used  in  the  tests  (fig.  2(c)).  A 
calibrated  five-tube  survey  rake  was  used  to  measure  the  stream  pitch 
and  yaw  angles  and  dynamic  pressures  for  most  of  the  tests  (fig.  2(d)). 
This  rake  is  fully  described  in  reference  12.  In  order  to  obtain  data 
closer  to  the  plane  of  the  rotor,  a  calibrated  pitch  head  (fig.  2(e)) 
was  used  to  survey  points  in  the  longitudinal  plane  of  symmetry  for  one 
flight  condition. 

The  procedure  used  in  the  tests  was  as  follows:  The  rotor  was  first 
set  approximately  at  a  predetermined  flight  condition.  Then,  the  area 
above,  below,  and  behind  the  rotor  was  surveyed  with  the  rakes.  At 
intervals  during  the  surveying,  readings  were  taken  of  rotor  thrust  and 
drag,  blade  flapping  and  feathering  motions,  and  tunnel  dynamic  pressure. 


CORRECTIONS 


A  floor  was  installed  in  the  tunnel  test  section  for  all  the  tests, 
and  for  this  configuration  the  Langley  full-scale  tunnel  has  no  appreciable 
jet-boundary  correction  (ref.  13). 

In  order  to  correct  the  measured  dynamic  pressures  and  stream  angles 
for  the  flow  about  the  supporting  mechanism  and  for  tunnel  stream  angle, 
data  were  taken  twice  at  each  point  in  space  —  once  with  and  once  without 
the  rotor  blades  installed.  The  difference  is  considered  to  be  charac¬ 
teristic  of  the  isolated  rotor.  It  should  be  noted,  however,  that  this 
procedure  is  not  adequate  in  regions  immediately  behind  the  rotor.  (See 
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figs.  3  "to  5*)  Measurements  made  in  these  regions  have  little  accuracy 
"because  of  the  effect  of  the  rotor  on  the  flow  about  the  tower.  The 
angle  of  attack  of  the  rotor  was  corrected  for  tunnel  stream  angle. 


PRECISION 


The  accuracy  of  the  two  rakes  is  essentially  the  same.  Both  measure 
dynamic  pressure  within  2  percent  and  stream  angles  within  l/2°.  These 
figures  take  into  account  the  unsteadiness  of  the  flow  and  the  high 
vibration  level  of  the  rake  mount.  At  the  lowest  tip-speed  ratio,  addi¬ 
tional  inaccuracy  is  caused  by  the  difficulty  in  reading  the  small  changes 
of  liquid  height  in  the  manometer  tubes.  Overall  accuracy  in  terms  of 
the  average  value  of  the  normal  component  of  induced  velocity  will  be  a 
function  of  the  rotor  lift  coefficient.  For  ease  of  reference,  the  fol¬ 
lowing  table  gives  the  estimated  overall  accuracy  for  each  of  the  flight 
conditions  covered  in  this  investigation: 


ft 


Tip-speed  ratio,  n 

Accuracy,  percent  Vq 

0.095 

±15 

.139 

±10 

.iko 

±10 

.232 

±25 

The  rakes  used  were  incapable  of  measuring  instantaneous  variations 
of  induced  velocity  such  as  might  be  desired  for  a  refined  blade -vibration 
analysis.  Therefore,  all  data  are  time-averaged  values.  The  experimental 
results  of  reference  8  indicate  that  instantaneous  values  may  differ 
greatly  from  the  mean  values  measured  in  these  tests. 

During  the  tests,  the  rotor  tip  speed  was  held  constant  to  within 
1  percent.  The  thrust  coefficient  was  measured  within  1  percent.  The 
skew  angle,  which  was  computed  by  using  the  measured  values  of  and 

the  measured  blade  motions  in  equation  (29)  of  reference  2,  should  be 
accurate  to  within  1°. 

The  position  of  the  rake  with  respect  to  the  hub  was  known  within 
0.3  percent  R.  Because  the  rotor  was  of  the  teetering  type,  the  only 
coning  was  due  to  the  deflection  of  the  blades.  The  tip  deflection  was 
estimated  visually  to  be  between  2  and  3  percent  of  the  radius.  The 
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data  are  presented  with  reference  to  a  plane  parallel  to  the  tip-path 
plane  and  passing  through  the  teetering  pin. 


RESULTS  AND  DISCUSSION 
Basic  Data 


Because  of  the  present  lack  of  information,  both  theoretical  and 
experimental,  on  the  induced  velocities  in  many  areas  near  the  rotor, 
substantially  all  the  basic  data  obtained  in  the  investigation  are  pre¬ 
sented.  In  each  case  the  rotor  flight  condition  is  identified  on  the 
figure  by  the  skew  angle  X  and  the  tip-speed  ratio  [i.  Table  I  iden¬ 
tifies  the  various  flight  conditions  more  completely. 

Figures  3  to  5  show  the  measured  stream  angles  in  the  flow;  fig¬ 
ures  6  to  8  show  contours  of  dynamic -pres sure  ratio;  and  figures  9  to  29 
give  the  measured  values  of  induced-velocity  ratio.  The  fairings  of 
the  measured  data,  shown  in  figures  9  to  29,  are  used  for  the  subsequent 
analysis . 

The  basic-data  curves  are  not  discussed  as  such  in  this  paper  since 
the  more  significant  phenomena  have  already  been  pointed  out  in  refer¬ 
ence  10. 


Discussion  of  Applicable  Theory 

In  reference  2,  Castles  and  De  Leeuw  present  calculated  values  of 
the  normal  component  of  induced  velocity  in  the  longitudinal  plane  of 
symmetry  and  along  the  principal  (X,Y,Z)  axes  of  a  rotor  with  a  uniform 
disk  loading.  Reference  3  uses  the  flow  field  of  reference  2  to  calcu¬ 
late  the  flow  field  for  rotors  having  nonuniform  but  circularly  symmet¬ 
rical  disk  loadings  (the  specific  cases  calculated  being  rotors  with  a 
triangular  disk  loading  and  with  a  typical  measured  mean  disk  loading). 

Additional  theoretical  calculations  by  Mangier  and  Squire  (ref.  4) 
give  induced-velocity  information  in  the  entire  rotor  disk  and  also  for 
a  transverse  plane  in  the  far  wake.  A  later  section  discusses  a  method 
of  transforming  these  results  to  correspond  with  the  results  of  refer¬ 
ences  2  and  3  except  for  the  specific  disk  loadings  considered. 

Drees  (ref.  5)  discusses  some  of  the  effects  of  a  finite  tip-speed 
ratio  and  presents  calculations  of  some  of  these  effects  at  selected 
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points  on  the  rotor  disk.  It  should  be  noted  that  the  calculated  curves 
presented  with  the  data  are  not  modified  for  the  asymmetries  discussed 
in  this  reference. 


Comparison  of  Measured  Data  With  Calculated 
Values  of  References  2  and  3 

Longitudinal  plane  of  symmetry.-  Figures  30  to  33  present  the  meas¬ 
ured  values  of  induced  velocity  in  the  longitudinal,  plane  of  symmetry  as 
compared  with  the  calculated  values  from  references  2  and  3« 

For  the  four  flight  conditions,  the  measured  data  near  the  forward 
portions  of  the  rotor  fall  close  to  the  flow  calculated  for  the  rotors 
with  zero  load  in  the  center.  The  existence  of  essentially  zero  induced 
velocity  at  the  center  of  the  rotor  and  of  regions  of  upwash  behind  the 
center  are  clearly  shown.  Further  rearward  in  the  flow,  just  behind  the 
trailing  edge  of  the  rotor  (x/R  =  1.07),  there  is  no  longer  any  pronounced 
correlation  between  the  measured  and  the  calculated  flow,  regardless  of 
the  lift  distribution  assumed  to  exist  on  the  rotor  disk.  At  greater 
distances  behind  the  rotor  (X/R  =  2.07  and  3*1^),  there  seems  to  be  no 
correlation.  The  physical  reasons  for  this  lack  of  correlation  are 
discussed  subsequently. 

Longitudinal  axis.-  If  curves  are  faired  through  the  data  in  fig¬ 
ures  30  to  33  and  the  values  at  Z/R  =  0  are  taken  from  these  curves, 
it  is  possible  to  obtain  an  idea  of  the  variation  of  induced  velocity 
along  the  longitudinal  axis  of  the  rotor.  The  results,  plotted  in  fig¬ 
ure  jb,  show  clearly  the  good  correlation  of  the  nonuniformly  loaded 
rotor  calculations  in  the  flow  over  the  forward  three-quarters  of  the 
rotor  and  the  deviation  behind  that  point. 

Figure  35  compares  the  measured  points  of  figures  34(c)  and  34(d). 
These  two  flight  conditions  have  very  small  differences  in  skew  angles. 
There  is,  however,  an  appreciable  difference  in  tip-speed  ratio  and 
lift  coefficient.  Actually,  the  calculated  induced-velocity  ratios  for 
these  two  skew  angles  are  virtually  identical.  It  may  be  seen  that, 
experimentally,  there  is  very  little  difference  between  the  induced 
velocities  along  the  longitudinal  axis  for  the  two  tip-speed  ratios. 

This  fact  is  as  assumed  in  references  2  and  5* 

Lateral  axis.-  The  basic  data  may  be  cross-plotted  in  the  same 
manner  to  obtain  the  induced-velocity  distribution  along  the  lateral  (Y) 
axis.  The  results  of  these  cross  plots  are  shown  in  figure  36  and  are 
compared  with  the  calculated  distributions  of  references  2  and  3*  The 
calculated  flow  is  symmetrical  since  no  account  is  taken  of  dissymmetries 
due  to  the  finite  tip-speed  ratio. 
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A  word  of  caution  is  necessary  at  this  point.  Because  of  the  neces¬ 
sity  of  maintaining  a  clearance  between  the  survey  rake  and  the  rotor, 
it  was  not  possible  to  obtain  data  much  closer  than  20  percent  R  above 
the  rotor.  (See  figs.  30,  32,  and  33-)  The  theory  predicts  rather  large 
•gradients  of  induced  velocity  at  locations  closer  to  the  rotor  than  this 
point.  In  the  longitudinal  plane  of  symmetry  it  is  possible  to  cross¬ 
plot  the  data  with  some  confidence  by  comparing  the  trends  shown  well 
above  the  rotor  by  the  theory  and  by  the  measured  points.  At  other  loca¬ 
tions,  without  the  theory  as  a  guide,  the  accuracy  of  cross  plots  is 
problematical.  Therefore,  the  following  discussion  of  the  induced -velocity 
distribution  on  the  Y-axis  should  be  regarded  as  being  on  a  less  firm 
basis  than  the  discussion  in  the  preceding  sections. 

Figure  36  indicates  that  the  rotor  does  have  zero  induced  velocity 
at  its  center.  This  is  in  agreement  with  the  assumption  that  the  rotor 
carries  zero  load  at  its  center  and  further  indicates  that  this  fact 
must  be  used  if  the  flow  is  to  be  calculated  with  any  degree  of  accuracy. 

Figure  36  also  shows  that  the  induced- velocity  peak  on  the  retreating 
side  is  higher  and  further  outboard  than  the  peak  on  the  advancing  side. 
This  result  is  logical  when  the  difference  in  loading  on  the  two  sides 
of  the  rotor  are  compared.  The  resultant  velocities  at  the  inner  portion 
of  the  blades  on  the  retreating  side  are  very  low  or  even  negative,  whereas 
the  corresponding  velocities  on  the  advancing  side  are  much  higher.  There¬ 
fore,  the  loading  on  the  inner  portion  of  the  retreating  side  of  the  disk 
must  be  less  than  that  on  the  advancing  side.  The  outer  portion  of  the 
retreating  side  must  then  carry  a  greater  loading  in  order  that  the  blade 
thrust  moment  remain  essentially  constant  around  the  disk.  Thus,  the 
loading  on  the  retreating  side  will  peak  farther  outboard  and  at  a  higher 
value  than  the  loading  on  the  advancing  side  of  the  disk.  This  is  the 
same  trend  shown  by  the  measured  induced  velocities  and  is  indicative  of 
the  very  close  relation  between  the  induced  velocities  and  the  loading 
along  this  axis. 

Comparison  of  Measured  Data  With  Calculated  Values  of  Reference  4 

Modification  of  calculations  for  finite  lift  coefficient.-  Mangier 
and  Squire,  in  reference  4,  compute  the  induced  velocities  in  the  rotor 
disk  and  in  a  transverse  plane  in  the  far  wake.  Their  calculated  results 
are  limited  to  vanishingly  small  lift  coefficients  in  that  the  trailing- 
vortex  system  leaves  the  rotor  in  the  same  direction  and  with  the  same 
velocity  as  the  free  stream. 

The  assumption  of  a  finite  lift  coefficient  alters  this  picture  in 
two  ways:  First,  the  trailing-vortex  system  must  now  cross  the  free- 
stream  flow  and,  second,  there  will  be  a  change  in  the  vortex  density. 
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The  results  of  reference  4  are  easily  modified  to  account  for  these 
changes  and,  when  so  modified,  will  correspond  with  the  results  of  ref¬ 
erences  2  and  3  except  for  the  specific  disk  loadings  considered.  The 
first  modification  is  merely  a  rotation  of  the  field  to  account  for  the 
inclination  of  the  wake.  The  incidence  angle  i  is  then  the  complement 
of  the  skew  angle  X.  The  change  in  vortex  density,  although  affecting 
the  absolute  induced  velocity  v,  will  have  no  effect  upon  the  induced- 


velocity  ratio  =  4 


0 


'(^l) 


Therefore,  the  new  value  of  v  may  be 


found  by  using  the  value  of  v/vq  from  reference  4  and  by  calculating 
from  an  equation  which  considers  the  finite  lift  coefficient  (such 
eq.  (32)  of  ref.  2). 


v0 

as 


For  the  purpose  of  comparison  with  the  measured  data  of  this  paper, 
all  the  induced  velocities  of  reference  4  have  been  modified  to  corre¬ 
spond  with  reference  2  in  regard  to  both  skew  angle  and  vortex  density. 


Plane  of  the  rotor.-  Values  of  induced- velocity  ratio  in  the  plane 
of  the  rotor  obtained  by  cross -plotting  the  measured  data  are  compared 
with  the  corrected  calculated  values  from  reference  4  in  figures  37  to  43. 
These  figures  are  discussed  only  in  general,  since  they  indicate  the  same 
results  as  the  previous  comparison  with  references  2  and  3-  This  result 
might  be  expected  since  the  calculated  results  should- be  comparable  if 
the  same  load  distributions  are  assumed.  The  resemblance  between  the 
typical  loading  of  reference  3  and  pressure  distribution  III  of  refer¬ 
ence  4  is  evident. 

In  general,  the  same  extent  of  agreement  is  noticeable  (see  especially 
fig-  37(b));  that  is,  the  calculated  flow  field  is  close  to  the  measured 
field  only  over  the  forward  three-quarters  of  the  disk. 

Each  figure  (figs.  37  to  43)  shows  that  the  induced- velocity  field 
for  pressure  distribution  III  (labeled  P.D.  Ill)  is  closer  to  the  measured 
field  on  both  sides  of  the  disk  than  that  for  pressure  distribution  I 
(labeled  P.D.  i).  Since  pressure  distribution  III  represents  a  disk 
loading  with  zero  load  in  the  center  and  the  other  distribution  does  not, 
this  result  is  as  expected. 


This  result  does  not  agree  with  that  of  reference  9j  which  indicates 
that  pressure  distribution  I  should  be  used  for  calculating  the  induced 
velocities  over  the  advancing  side  and  pressure  distribution  III,  over 
the  retreating  side  of  the  rotor.  This  result  (ref.  9)  arose  from  meas¬ 
urements  made  at  X/R  =  1.5  and,  as  such,  is  subject  to  error  because 
of  the  discrepancy  between  theory  and  the  measured  data  in  this  region. 
The  present  result  is  considered  to  be  more  accurate. 

The  values  at  X/R  =0.5  (fig-  ^3)  seem  to  indicate  that,  at  this 
location,  the  error  is  generally  somewhat  larger  near  the  outer  edges 


1 


I 


J 


\ 


10 


NACA  TN  3691 


of  the  disk.  This  result  might  be  expected  since  the  roll-up  of  the 
trailing-vortex  system  would  be  expected  to  affect  the  outer  edges  of  the 
disk  sooner  than  the  center. 

Far  wake.-  Reference  4  also  computes  the  induced-velocity  distribu¬ 
tion  in  the  far  wake.  The  computed  contours  of  induced-velocity  ratio 
for  both  pressure  distributions  in  the  far  wake  are  compared  in  fig- 
vires  44  to  49  with  measured  contours  of  induced  velocity  at  X/R  =  2.07 
and  X/R  =  3- l4  for  three  different  flight  conditions.  (Notice  that 
the  rotation  of  the  calculated  flow  field  required  to  obtain  correspond¬ 
ence  with  references  2  and  3  results  in  a  vertical  displacement  of  the 
computed  induced  velocities.) 

In  general,  the  measured  induced  velocities  are  of  the  same  order 
of  magnitude  as  the  calculated  values.  This  result  would  be  expected, 
regardless  of  the  flow  pattern  assumed  for  the  calculations,  since  a 
given  amount  of  downward  momentum  must  be  imparted  to  the  air  to  obtain 
a  given  lift.  The  distribution  of  induced  velocity  is,  however,  rather 
different  from  the  calculated  values,  particularly  in  the  central  por¬ 
tion  of  the  wake.  In  some  cases,  the  calculated  induced-velocity  field 
in  the  usual  location  of  short-span  tail  surfaces  could  be  incorrect  by 
as  much  as  1.5v/vq. 


Nature  of  the  Flow 

At  this  point  the  reason  for  the  rather  severe  discrepancies  between 
the  calculated  and  the  measured  induced  velocities  in  the  rearward  por¬ 
tion  of  the  flow  is  considered.  Reference  10  has  already  noted  the 
rolling  up  of  the  vortex  sheet  behind  the  rotor.  Figures  50  to  52  examine 
this  phenomenon  more  thoroughly  for  three  different  flight  conditions. 

These  figures  are  contour  maps  of  the  measured  local  mean  vorticity 
behind  the  rotor.  The  intersection  of  the  outermost  edge  of  the  assumed 
wake  (refs.  2,  3*  and  6)  with  the  survey  plane  is  shown  in  each  case 
by  the  dashed-line  ellipse.  If  the  assumed  wake  vortex  pattern  is  correct, 
the  contours  of  vorticity  should  be  along  this  wake  edge  (or  within  it, 
depending  upon  the  disk  load  distribution). 

In  general,  this  is  not  true  since  the  dominant  feature  in  each  case 
is  the  vorticity  associated  with  the  well-developed,  rolled-up  vortices 
behind  the  rotor  tips.  It  should  be  noted  that  these  vortices  have  only 
been  displaced  vertically  about  half  as  far  as  the  center  of  the  wake. 

The  lower  half  of  the  elliptical  sheet  is  essentially  completely  rolled 
up  at  the  trailing  edge  (figs.  50,  5l(a)>  52).  The  upper  half  of  the 

elliptical  sheet  (originating  from  the  rear  of  the  disk)  is  still  visible 
at  X/R  =  1.07>  but  it  also  is  completely  rolled  up  shortly  after  that 
point  (figs.  51(b)  and  51(c)). 
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The  similarity  between  the  flow  behind  the  rotor  and  that  behind  a 
low-aspect-ratio  wing  is  very  close  and  it  is  possible  to  draw  a  somewhat 
more  quantitative  analogy. 

Reference  l4  examines  the  rolling  up  of  the  trailing- vortex  system 
behind  a  low- aspect -ratio  wing.  For  lift  coefficients  of  the  order  of 
several  tenths,  the  roll-up  is  essentially  complete  in  a  very  short  dis¬ 
tance.  Figure  9  of  reference  l4  shows  that,  for  an  untwisted,  rectangular 
wing  of  aspect  ratio  4/n  (the  same  aspect  ratio  as  a  rotor),  and  for 
the  lift  coefficients  used  herein  ( 0.720,  0.373>  and  0.122),  the  roll-up 
distances  should  be  0.33>  0.64,  and  1-97  times  the  chord  length,  respec¬ 
tively.  It  is  more  difficult  to  determine  the  rate  of  roll-up  for  some¬ 
thing  as  complex  as  a  rotor,  but  figures  50  to  52  indicate  qualitatively 
that  it  is  also  extremely  rapid. 

These  observations  are  further  substantiated  by  the  induced- velocity 
measurements  previously  examined  in  figures  30  and  32.  The  measured  data 
in  these  figures  show  a  marked  discontinuity  (which  would  be  expected  at 
the  edge  of  the  wake)  only  in  the  vicinity  of  the  upper  wake  edge  at 
X/R  =  1.07.  Farther  behind  the  rotor  (X/R  =  2.07  and  3*l4),  there  is  no 
pronounced  discontinuity;  this  condition  indicates  the  absence  of  a  strong 
vortex  sheet  in  this  region. 

Figures  33(f)  and  33(g)  seem  to  show  a  discontinuity  somewhat  above 
the  upper  edge  of  the  assumed  wake.  The  previous  analogy  with  low-aspect- 
ratio  wings  indicates  the  possibility  that  the  roll-up  would  not  be  as 
rapid  in  this  case  because  of  the  low  rotor  lift  coefficient  (Cl  =  0.122). 

The  difference  between  the  calculated  and  the  measured  flow  fields  may 
be  explained  on  this  basis.  In  the  forward  portions  of  the  flow,  the 
elliptical  wake  has  not  rolled  up,  and  it  therefore  resembles  the  wake- 
vortex  model  used  in  the  calculations.  The  wake  here,  has,  of  course, 
the  greatest  effect  on  the  induced  velocities  in  the  forward  portions 
of  the  flow  and  essentially  determines  their  magnitude.  Consequently, 
in  this  part  of  the  flow,  calculations  give  a  reasonable  picture  of  the 
actual  induced  velocity. 

Farther  rearward,  however,  the  distorted  portion  of  the  wake  has  a 
greater  effect  in  determining  the  actual  induced  velocity.  The  calculated 
flow  field,  therefore,  is  increasingly  more  inaccurate  as  its  distance 
downstream  increases. 

A  short  distance  behind  the  rotor,  it  should  be  possible  to  assume 
that  the  flow  will  be  the  same  as  that  for  an  equivalent  wing.  The  min¬ 
imum  distance  behind  the  rotor  at  which  this  procedure  will  be  valid  is 
.a  function  of  the  lift  coefficient. 

From  reference  15  the  induced  velocity  in  the  plane  of  symmetry  of 
a  uniformly  loaded  rectangular  wing  is 
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The  flow  calculated  by  this  equation  is  symmetrical  about  Z/R  =  0; 
that  is,  it  does  not  allow  for  the  vertical  displacement  of  the  wake. 

As  pointed  out  previously,  the  centers  of  the  vortices  are  displaced 
downward  an  amount 


A  R  =  2  R  cot  x  (7) 

The  induced  velocity  for  this  assumed  wing  has  been  calculated  from 
equations  (6)  and  (7)  and  is  compared  with  the  flow  measured  at  X/R  =  2.07 
and  3.14  (for  p.  =  0.l4o  and  x  =  82.3°)  in  figure  53*  As  expected  in 
these  rearward  planes,  it  is  a  closer  approximation  to  the  actual  flow 
than  are  the  flows  calculated  on  the  basis  of  a  skewed  cylindrical  wake 
(fig.  32,  for  example). 

Another  point  is  worth  noting  here.  It  has,  upon  occasion,  been 
assumed  that  the  rotor  wake  must  have  more  vorticity  behind  the  retreating 
side  of  the  disk  than  behind  the  advancing  side  of  the  disk,  since  the 
unequal  blade  velocities  on  the  two  sides  require  the  retreating  blade 
to  have  a  greater  value  of  circulation.  This  would  be  true  only  at  the 
trailing  edge  of  the  blade  to  an  observer  riding  around  with  it.  However, 
to  an  observer  immediately  behind  the  rotor  disk  (as  in  these  surveys), 
the  summation  of  vorticity  behind  the  disk  must  be  equal  on  each  side, 
since  vorticity  can  not  end  in  space.  More  complete  contours  of  local 
mean  vorticity  than  those  given  in  figures  50  to  52  indicate  that  this 
is  essentially  true. 


CONCLUSIONS 


The  results  of  this  wind-tunnel  investigation  of  the  induced  flow 
near  a  single  rotor  in  simulated  cruising  and  high-speed  flight  are  as 
follows : 

1.  As  far  rearward  as  three-quarters  of  a  diameter  behind  the 
leading  edge  of  the  rotor  disk,  the  normal  component  of  induced  velocity 
of  a  lifting  rotor  may  be  calculated  with  good  accuracy  by  available 
theory  provided  that  a  realistic  nonuniform  disk  loading  is  assumed. 

2.  Rearward  of  this  three-quarter  diameter  point,  calculation  of  the 
induced  velocity  becomes  increasingly  inaccurate  because  the  trailing- 
vortex  system  rolls  up  very  rapidly,  and  the  wake  vortex  pattern  behind 
the  rotor  is  altered  from  that  assumed  as  a  basis  for  the  calculations. 

3.  At  locations  well  behind  the  rotor,  the  induced  flow  may  be  cal¬ 
culated  more  accurately  by  assuming  the  rotor  to  be  a  uniformly  loaded 
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rectangular  wing,  provided  that  the  rotor  lift  coefficient  is  sufficiently 
high. 


There  is  little  effect  of  the  as symmetries  associated  with  tip- 
speed  ratio  on  the  induced  flow  along  the  longitudinal  axis  of  the  rotor. 

5-  The  induced-velocity  calculations  of  Mangier  and  Squire  can  be 
modified  to  account  for  the  wake  crossing  the  flow  because  of  the  finite 
lift  coefficient. 


Langley  Aeronautical  Laboratory, 

National  Advisory  Committee  for  Aeronautics, 
Langley  Field,  Va.,  February  23,  1956. 
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Figure  2.-  Continued 
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Figure  2.-  Concluded. 
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vectors  represent  yaw  and  pitch  angles.  Base  of  vector  is  at  point  of  measurement.  X  =  75*0°; 

H  *  0.095. 
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Figure  3»-  Continued 
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Figure  3.-  Continued 
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Figure  3*-  Continued. 
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Figure  6.-  Contours  of  dynamic  pressure  q/^  in  the  vicinity  of  the  rotor 
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Figure  6.-  Continued 


Advancing  at  da 


Ratraattng  alia 


{tot ratting  ild* 


56 


NACA  TN  5691 


Advancing  aide  Retreating  side 


(a)  Z/R  *  0.52. 


(c)  Z/R  =  0.18. 


1.0  ,5  0  -.5  -l.o 

Y/R 

(d)  Z/R  =  0.11. 

Figure  9*-  Measured  values  of  induced-velocity  ratio  v/v  .  X/R  =  -1.0} 

X  =  75.0°}  n  =  0.095. 
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Figure  11.-  Concluded. 
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(f)  Z/R  =  -0.12. 
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Figure  12.-  Concluded. 
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Figure  13.-  Measured  values  of  induced-velocity  ratio  v/vQ.  X/R  =  1.07. 
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Figure  Ik.-  Measured  values  of  induced-velocity  ratio  v/vq.  X/R  *  2.07} 

X  =  75-0°}  n  =  0.095. 
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Figure  14.-  Continued. 
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Figure  15.-  Continued. 
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(g)  Z/R  =  -0.20. 


(h)  Z/R  =  -0.27. 
Figure  16.-  Concluded. 
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Figure  19. -  Concluded. 


NACA  TN  3691 


87 


KJi 

f 


3 


(b)  Z/R  =  0.17. 
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(d)  Z/R  ■  0.03. 

Figure  20.-  Measured  values  of  induced-velocity  ratio  v/vQ.  X/R  =  1.07; 

X  =  82.3°;  [i  =  0.140. 
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Figure  20.-  Continued. 
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Figure  20.-  Continued. 
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Figure  21.-  Continued. 
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Figure  21.-  Continued. 
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Figure  22.-  Continued. 
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Figure  22.-  Continued. 
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Figure  22.-  Continued. 
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Figure  28.-  Concluded 


Figure  30.-  Comparison  of  measured  and  theoretical  values  of  Induced- 
velocity  ratio  v/vq  in  longitudinal  plane  of  symmetry.  X  =  75*  0° 


Figure  JO.-  Concluded. 
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Figure  31.-  Continued. 
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Figure  32.-  Comparison  of  measured  and  theoretical  values  of  induced 
velocity  ratio  v/vq  in  longitudinal  plane  of  symmetry,  x  =  82. 
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Figure  56.-  Comparison  of  measured  and  theoretical  values  of  induced' 
velocity  ratio  v/vq  along  lateral  center  line. 
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(c)  X  =  83.9°;  n  =  0.21 
Figure  38*-  Concluded. 
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(a)  X  -  75.0°;  \i  =  0.095- 

Figure  39--  Comparison  of  measured  values  of  induced-velocity  ratio  v/vq 

at  Y/R  =  ±0-50  with  theory  of  reference  4  modified.  Pressure  distri 
butions  I  and  III  denoted  by  P.D.  I  and  P.D.  III. 


Figure  39*-  Concluded 
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(a)  x  =  75-0°;  u  =  0.095. 

Figure  40.-  Comparison  of  measured  values  of  induced-velocity  ratio  v/vq 

at  Y/R  ■  tO. 75  with  theory  of  reference  4  modified.  Pressure  distri¬ 
butions  I  and  III  are  denoted  by  P.D.  I  and  P.D.  III. 
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(a)  x  =  75.0°;  n  =  0.095. 

Figure  40. -  Comparison  of  measured  values  of  induced-velocity  ratio  v/v^ 

at  Y/R  =  ±0.75  with  theory  of  reference  4  modified.  Pressure  distri 
butions  I  and  III  are  denoted  by  P.D.  I  and  P.D.  III. 


(b)  x-  82.3°;  t*  *  0.140.’ 

Figure  40.-  Continued. 
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(c)  x  *  83. 9° $  n  ■  0.252. 
Figure  40.-  Concluded. 
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(a)  X  =  75- CT;  p.  =  0.095. 

■e  41.-  Comparison  of  measured  values  of  induced-velocity  ratio  v/v0  at  X/R  =  -0.5  with 
Leory  of  reference  4  modified.  Pressure  distributions  I  and  III  are  denoted  by  P.D.  I 
id  P.D.  III. 
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(a)  x-  75. 0°;  \i  =  0.095. 

Figure  43.-  Comparison  of  measured  values  of  induced -velocity  ratio  v/vq 

at  X/R  «  0.5  vith  theory  of  reference  4  modified.  Pressure  distri¬ 
butions  I  and  III  denoted  by  P.D.  I  and  P.D.  III. 
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(b)  Theoretical,  pressure  distribution  I. 


(c)  Theoretical,  pressure  distribution  III. 

Figure  44.-  Comparison  of  measured  values  of  induced-velocity  ratio  v/^ 
at  X/R  =  2.07  with  theory  of  reference  4  modified.  X  =  75. 0°; 
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(c)  Ulieoretical,  pressure  distribution  III. 

Figure  46.-  Comparison  of  measured  values  of  induced- velocity  ratio  v/ 
at  X/R  =  2.07  with  theory  of  reference  4  modified.  X  =  82. 3° j 
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(c)  Theoretical,  pressure  distribution  III. 

Figure  47.-  Comparison  of  measured  values  of  induced-velocity  ratio  v/ 

at  X/R  *  5*14  with  theory  of  reference  4  modified,  x  =  82.5°; 
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(a)  Measured  data. 


(b)  Theoretical,  pressure  distribution  I. 


(c)  Theoretical,  pressure  distribution  III. 

Figure  48.-  Comparison  of  measured  values  of  induced- velocity  ratio  v/’ 

at  X/R  =  2.07  with  theory  of  reference  4  modified.  X  =  83. 9° J 

n  =  0.232. 
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(a)  Measured  data. 


(b)  Theoretical,  pressure  distribution  I. 


(c)  Theoretical,  pressure  distribution  III. 

Figure  49.-  Comparison  of  measured  values  of  induced-velocity  ratio  v/vQ 
at  X/R  »  3.14  with  theory  of  reference  4  modified.  X=  83.9  i 
u  »  0.232. 
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Figure  50.-  Distribution  of  vorticity.  X/R  =  1.07;  X  =  75.0°;  pi  =  0.095. 
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